Charge ordering and hopping in a triangular array of quantum dots 
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We demonstrate a mapping between the problem of charge 
ordering in a triangular array of quantum dots and a frus- 
trated Ising spin model. Charge correlation in the low tem- 
perature state is characterized by an intrinsic height field or- 
der parameter. Different ground states are possible in the 
system, with a rich phase diagram. We show that electronic 
hopping transport is sensitive to the properties of the ground 
state, and describe the singularities of hopping conductivity 
at the freezing into an ordered state. 

Ordering and phase transitions in artificial structures 
such as Josephson junction arrays Q and arrays of quan- 
tum dots Q have a number of interesting properties. 
One attractive feature of these systems is the control on 
the Hamiltonian by the system design. Also, the exper- 
imental techniques available for probing magnetic flux 
or charge ordering, such as electrical transport measure- 
ments and scanning probes, are more diverse and flexi- 
ble than those conventionally used to study magnetic or 
structural ordering in solids. There has been a lot of the- 
oretical H and experimental Q studies of phase transi- 
tions and collective phenomena in Josephson arrays, and 
also some work on the quantum dot arrays 

There is an apparent similarity between the Joseph- 
son and the quantum dot arrays problems, because the 
former problem can be mapped using a duality transfor- 
mation on the problem of charges on a dual lattice. These 
charges interact via the D = 2 Coulomb logarithmic po- 
tential, and can exhibit a Kosterlitz-Thouless transition, 
as well as other interesting first and second order transi- 
tions ||. In the quantum dots arrays, charges are cou- 
pled via the D = 3 Coulomb 1/r potential, which can 
be also partially screened by ground electrodes or gates. 
In terms of the interaction range, the quantum dot ar- 
ray system is intermediate between the Josephson array 
problem and the lattice gas problems with short range in- 
teractions, such as the Ising model and its varieties. From 
that point of view, an outstanding question is what are 
the new physics aspects of this problem. 

A very interesting system fabricated recently || is 
based on nanocrystallite quantum dots that can be pro- 
duced with high reproducibility, of diameters ~ 15 — 
lOOAtunable during synthesis, with a narrow size dis- 
tribution (< 5% rms). These dots can be forced to 
assemble into ordered three-dimensional closely packed 
colloidal crystals ||, with the structure of stacked two- 
dimensional triangular lattices. Due to higher flexibility 
and structural control, these systems are expected to be 



good for studying effects inaccessible in the more tradi- 
tional self-assembled quantum dot arrays fabricated us- 
ing epitaxial growth techniques. In particular, the high 
charging energy of nanocrystallite dots, in the room tem- 
perature range, and the triangular lattice geometry of the 
dot arrays || are very interesting from the point of view 
of exploring novel kinds of charge ordering. 

Motivated by recent attempts to bring charge carriers 
into these structures using gates, and to measure their 
transport properties Q, in this article we study charge 
transport in the classical Coulomb plasma on a trian- 
gular array of quantum dots. We assume that the dots 
can be charged by an external gate and that conducting 
electrons or holes can tunnel between neighboring dots. 

For drawing a connection to the better studied spin 
systems, it is convenient to map this problem on the 
classical Ising antiferromagnet on a triangular lattice 
(AIAFM). Without loss of generality, we consider the 
case when the occupancy of the dots is either 1 or 0, 
and interpret these occupancies as the spin "up" and 
"down" states. In this language, the gate voltage is rep- 
resented by an external magnetic field coupled to the 
spins. Also, more generally, any spatially varying elec- 
trostatic field is mapped on a spatially varying magnetic 
field in the spin problem. For example, the electrostatic 
field due to, e.g., charged defects, corresponds to the 
spin problem with a random field. The electron hop- 
ping (tunneling) between the dots corresponds to spin 
exchange transport. The long-range Coulomb interaction 
between charges gives rise to a long-range spin-spin cou- 
pling, which leads to somewhat different physical prop- 
erties than the short-range exchange interaction conven- 
tional for the spin problems. 

Another difference between the charge and spin prob- 
lems is that, due to charge conservation, there is no ana- 
log of spin flips. This certainly has no effect on the 
equilibrium statistical mechanics, because the spin en- 
semble with fixed total spin is statistically equivalent to 
the grand canonical ensemble. However, this is known 
to be important in the dynamical problem Jl()| . The two 
corresponding types of dynamics are the spin conserv- 
ing Kawasaki dynamics p"l| ] and the spin non- conserving 
Glauber dynamics jL2|, respectively. In the simulation 
described below we use the Kawasaki dynamics, involv- 
ing spin exchange processes on neighboring sites. 

The mapping between the charge and spin problems 
is of interest because of the following. If the interac- 
tion between the dots were of a purely nearest neigh- 
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bor type, the problem could have been exactly mapped 
on the AIAFM problem, which is exactly solvable ||. 
The AIAFM problem is known to have an infinitely de- 
generate ground state with an intrinsic "solid-on-solid" 
structure described by the so-called "height field" (see 
below). The height field represents the correlations of 
occupancy of neighboring sites, and can be thought of 
in terms of an embedding of the structure into a three- 
dimensional space. The higher-dimensional representa- 
tions of correlations in solids have also been found useful 
in a variety of frustrated spin problems p3| . The or- 
dering of electrons in triangular arrays of quantum dots, 
because of Coulomb coupling giving rise to a repulsive 
nearest-neighbor interaction, must be similar to that of 
the AIAFM ground states. It represents, however, a new 
physical system in which the height field will be strongly 
coupled to electric currents, which can make electronic 
transport properties very interesting. 

The Hamiltonian of the electrons on the quantum dots 
is given by 7Y c har g c + Htunnci + W sp in, where 7i c har g o de- 
scribes Coulomb interaction between charges = 0, 1 on 
the dots and coupling to the background disorder poten- 
tial 4>{r) and to the gate potential V g : 

^charge = \Y1 V ( T v)<b<b + ^( V S + <K r 0)ft C 1 ) 

The position vectors run over a triangular lattice with 



the lattice constant 



ft, 



and Tij = ri 



r,. The interaction 



V accounts for screening by the gate: 



e\Ti 



{2df 



(2) 



Here e is the dielectric constant of the substrate^, and d 
is the distance to the gate plate. The single dot charging 
energy ^^(0) = e 2 /2C is assumed to be high enough to 
maintain no more than single occupancy. 

Electron tunneling between neighboring dots is de- 
scribed by Wtunnel- We assume that the tunneling is inco- 
herent, i.e., is assisted by some energy relaxation mech- 
anism, such as phonons. Below we consider stochastic 
dynamics in which charges can hop between neighboring 
dots with probabilities depending on the potentials of the 
dots and on temperature. Also, we consider only charge 
states and ignore all effects of electron spin described by 
7i sp in, such as exchange, spin ordering, etc. 

To lay out the framework for discussing ordering in 
the charge problem, let us review here the main results 



1 In the case of spatially varying e the interactions can be 
more complicated. For example, if the array of dots is placed 
over a semiconductor substrate, one has to replace e — » (e + 
l)/2 in the expression (^|). 




FIG. 1. Typical charge configuration obtained in a Monte 
Carlo simulation for the filling fraction n — 1/2. Pairing of the 
triangles is revealed by erasing all frustrated bonds connecting 
sites with equal occupancy. 



for the AIAFM problem using the charge problem lan- 
guage (and assuming nearest neighbor interactions) . The 
ground state of this system has a large degeneracy [||, 
which can be understood as follows. For each triangular 
plaquette at least one bond is frustrated, for any occu- 
pancy pattern. To minimize the energy of the nearest 
neighbor interactions, i.e., to reduce the number of frus- 
trated bonds, it is favorable to arrange charges so that 
the triangles combine in pairs in such a way that each 
pair of triangles shares a common frustrated bond. The 
pairing of triangles can be described by partitioning the 
structure into rhombuses. Given a charge configuration, 
the corresponding rhombuses pattern can be visualized 
by erazing all frustrated bonds, as shown in FigJ^. 

Conversely, given a configuration of rhombuses, one 
can reconstruct the arrangement of charges in a unique 
way, up to an overall sign change. The use of the rep- 
resentation involving rhombuses is that it leads to the 
notion of a height field. Any configuration of rhombuses 
can be thought of in terms of a projection of a faceted 
surface in a 3D cubic lattice along the (111) direction. 
This surface defines lifting of the 2D configuration in the 
3D cubic lattice space, i.e., an integer- valued height field. 
The number of different surfaces that project onto the do- 
main of area A is of the oder e w ^, where w is a constant 
equal to the entropy of the ground state manifold per 
plaquette. 

At a finite temperature, the charge ordering with re- 
spect to the height field, i.e., to the pairing of triangles, 
may have some defects (see Fig.||). An elementary de- 
fect is represented by an isolated unpaired triangle sur- 
rounded by rhombuses, i.e., by paired triangles. This 
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FIG. 2. Typical charge configuration for the filling frac- 
tion n — 1/2. In this case, the temperature is warmer than 
in Fig.[j], and there are several defects ("screw dislocations") 
present in the system. 



defect has a topological character similar to a screw dis- 
location, because it leads to an ambiguity in the height 
field. This ambiguity is readily seen in Figj^, where by 
continuing the height field around an unpaired triangle 
one finds a discrete change in it upon returning to the 
starting point. 

There is no finite temperature phase transition in 
the AIAFM problem because of the topological defects 
present at a finite concentration at any temperature. 
However, since the fugacity of a defect scales as e~ v ^ a ^ T , 
and the defect concentration scales with the fugacity, the 
system of defects is becoming very dilute as T — > 0. As a 
result, the T = state is ordered, and belongs to the de- 
generate manifold of ground states without defects, i.e., is 
characterized by a globally defined height variable. Also, 
since the distance between defects diverges exponentially 
as T — > 0, there is a large correlation length for the T = 
ordering even when T is finite. This situation is described 
in the literature on the AIAFM problem as a T — crit- 
ical point (see recent article JlTj and references therein). 

In this work, we study the charge system with the in- 
teraction (0) by a Monte Carlo (MC) simulation of elec- 
tron hopping in equilibrium, as well as in the presence of 
an external electric field. We find that, although many 
aspects of the AIAFM physics are robust, the long-range 
character of the interaction (||) makes the charge problem 
different from the AIAFM problem in several ways. 

The states undergoing the MC dynamics are all charge 
configurations with no more than single occupancy (qi — 
0, 1) on a square patch N x N of a triangular array 
(N = 12 in Figs.[l,||,||). Periodic boundary conditions 
are imposed by defining the energy (Q) using the N x N 



charge configuration extended periodically in the entire 
plane. Also, we allow charge hopping across the patch 
boundary, so that the charge disappearing on one side of 
the patch reappears on the opposite side, consistent with 
the periodicity condition. 

The stochastic MC dynamics is defined by letting elec- 
trons hop on unoccupied neighboring sites with proba- 
bilities given by Boltzmann weights: 



= e (*« 



-<S> 3 )/kT 



Wi. 



+ W^i = 1 , (3) 



where 



(4) 



To reach an equilibrium at a low temperature, we take 
the usual precautions by running the MC dynamics first 
at some high temperature, and then gradually decreasing 
the temperature to the desired value. 

All the work reported below was done on the system 
without disorder, (f>{vi) = 0. The distance to the gate 
which controls the range of the interaction (||) was chosen 
to be d = 2. 

The properties of the system depend on the charge 
filling fraction, i.e., on the mean occupancy n = 1i/N 2 , 
which is conserved in the MC dynamics. (The AIAFM in 
the absence of external magnetic field corresponds to n = 
1/2.) Wc find that at low temperature the equilibrium 
states with 1/3 < n < 2/3 are very well described by 
pairing of the triangles, as illustrated in Figs.[l],|3|. The 
defects with respect to this height field ordering have a 
very small concentration, if present at all. The short- 
range ordering in the charge problem turns out to be the 
same as in the AIAFM. Qualitatively, this similarity is 
explained by a relatively higher magnitude of the nearest 
neighbor coupling (|J) compared to the coupling at larger 
distances. 

Like in the AIAFM, we observe topological defects. 
They are present at warmer temperatures (see FigM and 
quickly freeze out at colder temperatures (see Figjl])^ The 
disappearance of the defects, because of their topological 
character, is possible only via annihilation of the opposite 
sign defects. An example of such a process can be seen 
in the lower right corner of Fig.||. The defects in the 
charge problem, besides carrying topological charge, can 
carry an electric charge. Two defects of opposite electric 
charge can be found at the top and on the left of Fig.||. 

Expectedly, besides the robust features, such as the 
height field and the topological defects, there are certain 
non-robust aspects of the problem. The two most inter- 
esting issues that we discuss here are related with lifting 
of the degeneracy of the AIAFM ground state, and with 
the instability of the T = critical point, both arising 
due to long-range coupling in the charge problem. 

As we already have mentioned, the ground state man- 
ifold of the AIAFM is fully degenerate only for purely 
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Charge density per lattice site 



FIG. 3. Typical Monte Carlo charge configuration for 
n — 1/3 + e. There are three excess charges in the system 
(e = 3/144) hopping over the honeycomb network of unoccu- 
pied sites in the commensurate y/S x y3 state. The excess 
carriers, dilute at e < 1, are moving nearly independently on 
the background of the frozen \/3 x \/3 state. 

nearest neighbor interactions. The degeneracy is lifted 
even by a weak non-nearest-neighbor interaction |]l4] , |l~5| . 
For instance, for our problem with the interaction (J2|) 
at n = 1/2, the favored ground states have the form 
of stripes, spaced by V3 in the lattice constraint units, 
which corresponds to electrons filling every other lattice 
row. Because there are three possible orientations of the 
stripes, the system cooled in the n=l/2 state freezes in 
a state characterized by domains of the three types. One 
of such domains of stripes can be seen in the upper left 
part of Fig.|l|. Upon long annealing, the domains some- 
what grow and occasionally coalesce. However, we were 
not able to determine whether the system always reaches 
a unique ground state, or remains in a polycrystallinc 
state of intertwining domains. 

This can be contrasted with the behavior at n = 1/3, 
where the ground state is a \/3 x y3 triangular lattice, 
corresponding to electrons filling one of the three sublat- 
tices of the triangular lattice (see Fig.||). The rotational 
symmetry of this state is the same as that of the un- 
derlying lattice. Because there is a (triple) translational 
degeneracy of the ground state, but no orientational de- 
generacy, the system always forms a perfect triangular 
structure upon cooling, without domains. For the den- 
sities close to 1/3, the ground state contains \n — 1/3] 
charged defects, vacancies for n < 1/3 and interstitials 
for ii > 1/3, formed on the background of the otherwise 
perfect -s/3 x structure. At n > 1/3, the interstitials 
are moving over the honeycomb network of empty sites, 
as illustrated in Fig.|[ 



FIG. 4. Conductivity versus electron density for several 
temperatures, given in the units of e 2 /ea. The density n varies 
from (uncharged) to 1 (fully charged), however, because of 
the n <-> 1 — n symmetry, only the interval < n < 1/2 is dis- 
played. Arrows mark the features corresponding to the phase 
transitions at n — 1/3, 1/2, 1/4. 

Before discussing the issue of finite T versus T = 
phase transitions, let us explain how we use the MC simu- 
lation to find electrical conductivity. It is straightforward 
to add an external electric field E to the MC algorithm 
|i"6f . For that, one can simply modify the expressions (||) 
for the hopping probabilities by adding the field poten- 
tial difference E • (r^ ) between the two sites and rj. 
In doing this, one has to respect the periodic boundary 
condition, which amounts to taking the shortest distance 
between the sites and r, on the torus. Then the charges 
are statistically biased to hop along E, which gives rise 
to a finite average current J. 

The current J is Ohmic at small E, and saturates 
at |E|a > min(U(a)n 1 / 2 , kT), where a(n) = n^ 1 / 2 is 
the inter-electron spacing. Even though only the Ohmic 
regime J ~ E is of a practical interest, it is useful to 
maintain the field E used in the MC simulation on the 
level of not less than few times below the nonlinearity off- 
set field, to minimize the effect of statistical fluctuations 
on the time averaging of the current J. 

The dependence of the electric current on the density 
n is shown in Fig.|] for several temperatures. Because 
of the electron-hole symmetry of the system, the con- 
ductivity problem is invariant under the transformation 
n <-> 1 — n, and thus only the interval < n < 1/2 can 
be considered. At high temperatures kT > e 2 /ea(n), the 
Coulomb interaction is not important, and one can eval- 
uate the current by considering electrons freely hopping 
over lattice sites subject only to the single occupancy 
constraint. The current J in this case goes as 
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FIG. 5. Conductivity versus temperature for density 
n — 1/3. The three curves are obtained for three different 
system sizes. 

J = -n(l - n) — , (5) 

where the factor n(l — n) in (^|) gives the probability that 
a particular link connects two sites of different occupancy, 
so that hopping along this link can occur, whereas the 
factor (E/fcT) in (||) is the hopping probability bias due 
to the weak field |E| <C kT. The constant factor 3/4 
in (JsJ) is given by the coordination number of the lattice 
(equal to 6) divided by 8. 

The parabolic n(l — n) dependence of the current (|5|) 
is clearly reproduced by the highest temperature curve 
in Fig.^. In Fig.^j the inverse temperature factor of the 
expression (||) is eliminated by rescaling the current by 
E/kT. The rescaled current J(kT/E) in Fig.|| is prac- 
tically temperature independent at small densities n, 
in agreement with the result d). Regarding the 1/kT 
rescaling, note that in a real system the frequency of 
electron attempts to hop is a function of temperature, be- 
cause hopping is assisted by energy relaxation processes, 
such as phonons. In this case, our MC result for the cur- 
rent J has to be multiplied by some function A(T), which 
will be of a power law form T a for the phonon absorption 
rate in the case of a phonon-assisted hopping. 

At the temperatures of the order and smaller than 
e 2 /ea(n), the current becomes suppressed due to charge 
correlations reducing the frequency of hopping. The on- 
set of this suppression takes place at temperatures of 
order of the interaction between neighboring electrons, 
kT ~ 0.2V(a)n 1 ' 2 . Note that the onset temperature is 
lower for smaller density, in agreement with the plots in 
Fig.| 

The electrical conductivity is a sensitive probe of freez- 
ing transitions, at which the system of interacting charges 



locks collectively into a particular ordered or disordered 
state, and the conductivity vanishes. Whether the freez- 
ing occurs as a finite T or a T — phase transition is 
related to the degree of degeneracy of the ground state. 
The situation appears to be very different at different 
densities n. For simple rational n — 1/3, 1/2, 2/3, 
and alike, characterized by the ground state unique up 
to discrete symmetries, there is a well defined freezing 
temperature. This is illustrated in Fig.||, where electri- 
cal conductivity of the n = 1/3 state is plotted versus 
temperature. The abrupt drop of the conductivity at 
kT w 0.09e 2 /ea, where a is the lattice constant, indi- 
cates a sharp freezing transition. To make sure that this 
is not a finite size effect, we show the conductivity curves 
for systems of three different sizes, 6x6, 12 x 12, and 
24 x 24. 

On the other hand, away from the specific densities 
with simple ground states, the freezing appears to be 
very gradual, and in the temperature range we explored 
in Fig^ there has been no evidence of a sharp transition. 
It may be that the ordering actually takes place at much 
smaller temperatures than the characteristic interaction, 
the situation not uncommon in frustrated systems. Also, 
it may be that at incommensurate densities the state re- 
mains disordered down to T = 0. From our observations, 
the latter seems to be a more likely scenario. We expect 
that upon cooling the charges freeze into a quasirandom 
state determined by the cooling history. In that case, 
this system represents an electronic glass that exists in 
the absence of external disorder. 

The nature of the ground state in this case is unclear. 
To list several options, it may be that the system forms 
a polycrystal consisting of intertwining domains, like it 
does at n = 1/2, or that the state represents a distorted 
incommensurate charge density wave, or that it is a gen- 
uine glass. Studying this would require enhancing the 
MC algorithm to make it capable of treating the slow 
dynamics of annealing at low temperatures. 

In conclusion, correlations of charges in the triangu- 
lar arrays are described by an intrinsic order parameter, 
the height field, similar to the AIAFM problem. The 
charge-spin mapping shows that various interesting phe- 
nomena arising in frustrated spin systems can be studied 
in charge systems, for which more powerful experimental 
techniques are available. The type of order in the ground 
state depends on the charge filling density. Electron hop- 
ping conductivity is sensitive to charge ordering and can 
be used as a probe of the nature of the ordered state. 
Conductivity couples to the height field, and thus one 
can expect novel effects in electronic transport proper- 
ties which have no analog in spin systems with the same 
geometry. 
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